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Abstract 

We solve the loop equations of the /3-ensemble model analogously to the 
solution found for the Hermitian matrices = 1. For (3 = 1, the solution was 
expressed using the algebraic spectral curve of equation y 2 = U (x). For arbitrary 
j3, the spectral curve converts into a Schrodinger equation ((hd) 2 — U(x))ip(x) = 
with h oc {yfi3~—l/y/~j3)/N. This paper is similar to the sister paper I, in particular, 
all the main ingredients specific for the algebraic solution of the problem remain 
the same, but here we present the second approach to finding a solution of 
loop equations using sectorwise definition of resolvents. Being technically more 
involved, it allows to define consistently the ,8-cycle structure of the obtained 
quantum algebraic curve (a D-module of the form y 2 — U (x) , where [y,x] = h) and 
to construct explicitly the correlation functions and the corresponding symplectic 
invariants J-h, or the terms of the free energy, in 1/N 2 -expansion at arbitrary K. 
The set of "flat" coordinates comprises the potential times % and the occupation 
numbers e a . We define and investigate the properties of the A- and ^-cycles, 
forms of 1st, 2nd and 3rd kind, and the Riemann bilinear identities. We use 
these identities to find explicitly the singular part of J-q that depends exclusively 
on e a . 
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1 Introduction 

In the contemporary mathematical physics, one can often meet the notion of quantum 
surfaces, which appears in many different aspects. Having no intension to describe 
all problems in which quantization of the very space-time coordinates takes place 
(which pertains mainly to string or brane models) we however stress that the main 
feature of most, if not all, these models is that the consideration is commonly restricted 
to simple geometries of sphere or torus. Observables in these theories are not the 
coordinates, which cease to commute with each other and satisfy some postulated 
quantum algebras, but objects related to representations of these algebras, because 
only these objects admit classical interpretation. In this paper, we propose a new 
approach to the description of these so-called "quantum surfaces," namely we begin 
with solutions of the standard one-dimensional Schrodinger equation with a polynomial 
potential and construct a higher genus quantum surface (which is the analogue of a 
classical hyperelliptic Riemann surface) for which we can define analogues of all the 
main notions of algebraic geometry. 

This paper is an "alternative version" of our paper jS] in which the notion of the 
quantum algebraic geometry was introduced and which we refer to as paper I in what 
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follows. In the both versions, the origin of quantum algebraic geometry is the same, 
the Schrodinger equation ((hd) 2 — U(x))ip(x) = 0. The principal difference is that in 
this, second version, we use the sectorwise definition of all the quantities starting from 
the one-point resolvents, i.e., we use different solutions of the Schrodinger equation 
to construct these resolvents in different Stokes sectors of the complex plane. This 
enables us to define in a rigorous way the integrations over A- and i3-cycles as well 
as to present a self-consistent procedure for constructing the correlation functions and 
the symplectic invariants. 

The correlation functions Wn h \x\, • • • , x n ) and the symplectic invariants Fh for any 
algebraic plane curve given by a polynomial equation 

£(%,y) = J2 £i >i xiyj = 

were defined in [HJE2]- The invariants J-"h{£) are defined in terms of algebraic geometry 
quantities related to the Riemann surface of equation £{x, y) = 0. On the matrix model 
side, these invariants are terms of the 1/N 2 - (the genus) expansion of the free energy 
calculated in jl] for the one-matrix model and in [6] for the two-matrix model. 

We introduce the notion of a "quantum curve" for which S(x,y) is a non- 
commutative polynomial of x and y: 

£(v,y) = J^jJcV » [y,x] = h. (1.1) 

i,3 

The notion of quantum curve is also known as D-modules, i.e., a quotient of the space 
of functions by Ker S(x, y), where y = hd/dx. 

Our construction is based on the functions ij)(x) such that 

£{x,hd x ) = (1.2) 

and we show that one can consistently define all the basic notion of algebraic geometry 
within this construction. Whereas some objects, like branch points, become obsolete, 
we can define cycles, forms, Bergman kernel, period matrix and the corresponding Abel 
maps as well as other objects in a consistent way. 

But, otherwise, it is striking to find that almost all relationships of classical algebraic 
geometry remain unchanged when h ^ 0, for instance, the Riemann bilinear identity, 
the modified Rauch variational formula, and the topological recursion defining the 
correlation functions and the symplectic invariants. 

The symplectic invariants Th were first introduced for the solution of loop equations 
arising in the 1-hermitian random matrix model [9j 0] . They were later generalized to 
other hermitian multi-matrix models [U IT3] . 
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The models that correspond to the quantum surface are the /3-ensembles classified 
by the exponent (3. The three Wigner ensembles (see [T7], and we changed (3 — >■ (3/2) 
correspond to = 1 (hermitian case), (3=1/2 (real symmetric case), (3 = 2 (real 
self-dual quaternion case), but we can easily define a /3-ensemble eigenvalue model for 
any real value of (3 as the iV-fold integral of the form 

J d\t ■ ■ ■ d\ N \A(X)\^e- N ^^ v ^ 

(A is the Vandermonde determinant). 

In [3J, the solution of [9 J was generalized to the /3-ensembles, but the solution was 
presented as a double half- infinite sum for (3 = 0(1) at large N, 

oo 

F=J2 N2 ~ 2h ~ k (VP- VvW^M- (1.3) 

h,k=0 

The coefficients Th,k were computed in [3J. 

In this paper, as in paper I, we assume that h = (vT^ — VVTD/^j so we perform an 
(infinite) resummation in the above formula; the free-energy expansion then acquires 
the standard form, 

oo 

J= N 2 - 2h F h (h). (1.4) 

h=0 

The J-fi,k s of [3] can be recovered by computing the semi-classical small ^-expansion 
of J-h{K)- We demonstrate that J^hifi) is the natural generalization of the symplectic 
invariants of [12] for a "quantum spectral curve" S(x,y) with [y,x] = H. 
We define also analogues of the multi-point resolvents 

IN N 

W n (xi, ...,x n ) = N~ n I Y^(xi - A,)- 1 ■ ■ ■ Y,(xn - Xj)- 

\ 3=1 3=1 

where we let angular brackets denote the averaging with the weight 
|A(A)| 2/3 e _Arv/ ^^ =1 v ( x i\ These resolvents in turn admit the l/iV 2 -expansion, 
W n (x\, . . . , x n ) = XlhLo N 2 ~ 2h ~ n Wn h \xi, . . . , x n ), and we calculate all the terms wi^ 
using the modified diagrammatic technique. 

The main tool applied for studying the /3-eigenvalue model is the loop equation 
method. We obtain loop equations from the invariance of an integral under the special 
change of variables. Loop equations for the /3-eigenvalue model can be found in [7], [TO] , 
and here we solve them order by order in 1/N 2 , at fixed h. 

Recently, models of this type got a new vim due to the conjecture by Alday, Gaiotto, 
and Tachikawa (AGT) [TJ relating Nekrasov's instanton function [18] to conformal 
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blocks of the Liouville theory; these conformal blocks in turn can be described by 
the matrix-like model (see [16], [5]); the relation to the Nekrasov's parameters is 
explicit: eie 2 ~ 1/N 2 and e 1 /e 2 ~ 0, so using the approach in this paper, we can 
construct nonperturbative solutions of Nekrasov's formulas in ex/e^- in this paper, we 
investigate only the case of polynomial potentials, the generalization to the realistic 
logarithmic potentials appearing in the AGT conjecture will follows. 

The structure of the paper is as follows: we collect the generalities on the Stokes 
phenomenon pertaining to solutions of the Schrodinger equation in Sec. We describe 
our quantum Riemann surface in Sec. E] where we introduce A- and i3-cycles, filling 
fractions <Tj, and the first-kind functions (analogues of holomorphic and Krichever- 
Whitham meromorphic differentials) as well as the system of flat coordinates and the 
Riemann period matrix. In Sec. HJ we introduce the recursion kernels and the second- 
and third-kind (bi-)differentials. In Sec. EJ we go beyond the leading approximation 
in 1/N 2 and construct correlation functions of all orders using the Feynman-like dia- 
grammatic technique. We reveal the origin of our recursion procedure in Sec. El where 
we develop in details the variations w.r.t. the set of flat coordinates; the summary 
is in Sec. [71 In the next two (completely new as compared to paper I) sections, we 
investigate the link to the /3-ensemble models (Sec. E]) and construct on the base of 
this analysis the free-energy terms (Sec. |9]). In the first three appendixes to the paper, 
we present proofs of the three main theorems of Sec. |5] concerning properties of the 
correlation functions whereas the fourth appendix contains the new formula expressing 
J-q through the filling fractions e,; in the matrix model approach, the singular term 
has the structure |ef log q whereas in the quantum geometry this term is proportional 
to J logr(?j), which is the first actual example of calculations in the case of quantum 
Riemann surfaces. 

2 Schrodinger equation and resolvents 

2.1 Solutions of the Schrodinger equation 

We begin with the Schrodinger equation 



with U(x) being a polynomial of even degree 2c? for which we define the polynomial 
"potential" V{x) of degree d + 1 to be 



H 2 ip"(x) = U(x)^(x) 



(2.1) 



d 




(2.2) 



fc=0 
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We also define the polynomial of degree d — 1, 

P(x) = _ _ ^ (2 .3) 

Eventually, we define: 

xP(x) . 
to = km — (2.4) 

In the matrix model language (see section [SJ, £i, . . . , t^+i are called the times associated 
to the potential V(x), to is the normalized total number of eigenvalues (particles), or 
the temperature, whereas the remaining coefficients of P are defined by introducing 
fixed "filling fractions" below. 

2.1.1 Stokes Sectors 

A function ip{x) that is a solution of the Schrodinger equation exhibits the Stokes 
phenomenon, i.e., although ip(x) is an entire function, its asymptotics are discon- 
tinuous near oo where it has an essential singularity. Let On = Arg(£^+i) be the 
argument of the leading coefficient of the potential V(x). We define the Stokes 
half-lines the asymptotic directions along which KeV(x) vanishes asymptotically, 
Lk = <yX / Arg(x) = —jfj + 7r j, together with the corresponding Stokes sectors: 



S k = \ Arg(x) G 



$0 k 7T On k -\- \ 

h IT h IT 

d+l d+V d+1 d+l 



(2.5) 



i.e., Sk is the sector between Lk-\ and L*k- 

Notice that in even sectors we have asymptotically ReV(x) > and in odd sectors 
we have ReV(x) < 0. 

2.1.2 The Stokes phenomenon. Decreasing solution 

From the study of the Schrodinger equation it is known that ip(x) is an entire function 
having a large x expansion in each sector Sk, 

M X ) rsj e ± ^ v{x) x Ck (A k + — + ...) (2.6) 

s k x 

and the sign ±, may jump discontinuously from one sector to another as well as the 
numbers Ak, Bk,Ck, ■ ■ ■ (and in general, all the coefficients of the series in i at infin- 
ity)B 

In every sector Sk there exists a unique solution that decreases exponentially along 
each direction inside the sector. We now separate solutions in the even and odd sectors 



1 The corresponding series is asymptotic, so we cannot continue it analytically to other sectors. 
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no zero concentration 



Figure 1: Example of the Stokes sector partition and structure of zeros for the 
Schrodinger equation solution if>(x) that decreases in the light-colored sector and in- 
creases in all other sectors (the degree of the potential V(x) is four). 



and consider the set {ip a (x)} of solutions each of which decreases in the corresponding 
even sector. We therefore introduce a sectorwise system of solutions to the Schrodinger 
equation. 

An important and useful result is the Stokes theorem, which claims that if the 
asymptotics of ip(x) is exponentially small in some sector, then the same asymptotic 
series expansion (12.61) is valid in the two adjacent sectors (and therefore if>(x) is expo- 
nentially large in those two sectors). 

In the general case, (i.e., for a generic potential U(x)), the solution if> a (x) de- 
creases only in the sector S a , and is exponentially large in all other sectors. But if the 
Schrodinger potential U(x) is non-generic, then there may exist several sectors in which 
if> a {x) is exponentially small (which means that if> ai (x) = ip a2 (x) for some «i 7^ a^)- 

In what follows, we mainly consider the general case, so in what follows we assume 
all the functions if) a to be different if not stating the opposite. 

The case studied in [11] was the most degenerate case in which one and the same 
solution if) is exponentially small in d + 1 sectors. 

2.1.3 Zeroes of ip 

Every if> a (x) is an entire function with an essential singularity at 00, and with isolated 
zeroes s\ a \ ip a (s\ ) = 0. The number of these zeros can be finite or infinite. In the 
latter case, zeroes may only accumulate near 00, and only along the Stokes half-lines 
Lj bordering the sectors (see fig J2. 1.1ft . This accumulation of zeroes along the half-line 
Lj occurs if and only if if) a {x) is exponentially large on both sides of the half-line. So, 
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no accumulation of zeros of the function ip a (x) occurs along the lines that border the 
a's sector, and this function can therefore have only a finite number of zeros inside 
the "bigger" sector when we join the a's sector with the adjacent parts of the two 
neighboring sectors. 

If U(x) is generic, then each of ip a {x) has an infinite number of zeroes, the zeroes 
accumulate at oo along all half-lines Lj with j ^ a, a — 1. 

In paper I, we define the genus of the Schrodinger equation to be related to the 
number of half-lines of zeros accumulation of a selected function ipo. However, this 
definition is scheme- dependent, and we can in principle obtain different genera for the 
very same function U(x). The clear understanding of this is still lacking; a possible 
explanation is that we actually deal with different sections of an ambient infinite-genus 
quantum surface. 

2.1.4 Sheets 

In sector S a we have the asymptotic behavior 



and the function ip a has the same asymptotic behavior in the two adjacent sectors. 

We consider an a's sheet of the quantum Riemann surface to be the union of these 
three sectors with possible analytic continuation into a finite domain of the complex 
plane. We consider only the sheets enumerated by even a and, in contrast to paper I, 
introduce democracy of sheets: all of them will be equivalent in the approach of this 
paper. Sheets obviously overlap; we have to choose boundaries (cuts) between them. 

2.2 Resolvent 

The first ingredient of our strategy is to define a resolvent similar to the one in matrix 
models. 

Definition 2.1 We define the resolvent sectorwise: 



For a quantity defined sectorwise we indicate it by setting the sector index above 
the variable, as shown in (12. 8j) . 

It follows from this definition that u(x) has simple poles at zeros of ip a in the 
corresponding sector. The boundaries between sectors overlap, but in what follows we 
fix them in a more explicit form (see the partition of the complex plane by ^4-cycles). 




(2.7) 
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A straightforward computation then gives 

w(S) ~ - + 0(l/x 2 ), (2.9) 

that is, in each sheet the resolvent possesses asymptotic properties of a standard matrix- 
model resolvent. 

An important property of any solution i\) a is that 

Res— ^- = (2.10) 

The main property of u;(x) is that it satisfies the Ricatti equation. We obtain 

V"(x)u;(a;) - oj\x) - hu'(x) = — — fir ) { - fi — — 

4 ip a (x) 2 

= P(x), ' (2.11) 

with P(x) being a polynomial of degree d — 1 in x and this polynomial is one and the 
same for all sheets of the quantum Riemann surface introduced below. 



3 Quantum Riemann Surface 

In this section we define the notions of A- and £>-cycles and the first kind differentials 
dual to them. 

3.1 The contour Co and the set of A- and ^-cycles 

In papers on matrix models (on an early stage, before coming to residues at the branch 
points), we have the special contour of integration, Cd, that encircles all the singularities 
of resolvents leaving apart all other possible singular points. The analogue of such a 
contour in our case is the union of d + 1 contours, one per each sheet, that pairwise 
coincide in far asymptotic domains of odd Stokes sectors and separate all the zeros of 
the function ip a from the infinity oo a (which is always possible because we have a finite 
number of zeros in each sheet). We have (see Fig. 13. II) 

f f{x)dx = / f&dx (3.1) 
for any function fix) that has no asymptotic zero concentration along the boundary 

ex. 

lines of the sector S a . Here and hereafter, we assume that f[x) may depend on a finite 
number of derivatives of the function ip a (x), the symbol f[x) then indicates that we 
substitute the solution if) a (x) as an argument. 
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3.1.1 A- and i3-cycles 

We now deform the integration contour Co pushing through the "middle" part of the 
complex plane and taking the residues at the zeros of the corresponding functions 
ip a as shown in Fig. 13.1.11 On the way we might break some contours presenting them 
as the unions of newly introduced contours all of which are stretched between different 
asymptotic directions. As a result, we obtain a system of exactly 2d contours in which 
(leaving aside the residues at zeros s^) all the contours are pairwise identified and 
represent edges of d "cuts". As the result, we obtain a complete system of d cuts 
Ai, i = 1, . . . , d, that separate all the odd- numbered infinitiecl and determining the 
corresponding sheets of the quantum Riemann surface. If the functions ip a ( x ) coincide 
for some sheets, then we can identify these sheets. Note that we definitely have an 
arbitrariness in constructing this system of cuts; we can also arbitrarily assign the 
residues inside the sheet to belong to one of several contours bounding this sheet. 

We call the cut separating two sheets a cycle A a , and it is characterized by four 
indices: a + and a_ are indices of the sheets separated by this cut (they are even 
numbered in our classification); a + and 5L are indices of infinities that are asymptotic 
for this cut (they are odd numbered). 

To each complete set {^4 Q }^, =1 of ^4-cycles we uniquely set into the correspondence 
the set {B a }a=i of i3-cycles that go pairwise between the even-numbered infinities (oo a+ 

2 In what follows, we identify an infinity "point" with the corresponding number with the related 
asymptotic direction. 
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and oOq, ) such that the intersection index A a o J3p — o~a,p- 

Definition 3.1 We define the integrals over the cycles A a and the conjugate cycle B a 
to be (see Fig. 13X3}) 

/ f(x)dx= [°° a+ (fCx)-fCx))dx + J2 res f(°x) (3.2) 

and 

f f(x)dx= / (fCx)-f( a £))dx, (3.3) 

J B a Joo a _ 

where the residues in the first expression are taken at those zeros of ip a± that are 
assigned to the corresponding contour. 

Because the prescription for the sheet assignment follows from the definitions (13.21) 
and (13.31) of the cycle integrals, we omit the sheet labels in the corresponding integrands. 




Figure 3: Example of pushing the contour Cd from infinities to the set of ^4-cycles. 



Remark 3.1 Assignment of residues in the ath sheet to the contours bounding this sheet 
is arbitrary; we have therefore a (discrete) ambiguity in the definition (|3.2j) of the ^4-cycle 
integrals. However, the notion of the integral over Cd is well defined and does not depend on 
the choice of ^.-cycles. 

Obviously, § Cd f(x)dx = £)f =1 f A . f(x)dx. 

We now introduce the "genuine" A- and £>-cycles, which are straightforward ana- 
logues of the set of A- and i3-cycles on a standard Riemann surface. For this, we 
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select one among the ^4-cycles, say, the cycle Ad and the conjugate cycle Bd- Then, we 
identify A% = Ai and B{ — Bi — Bd for i — 1, . . . , d — 1 in the sense of Definition 13.11 
that is 

(p f(x)dx = f <p f(x)dx, 

(p f(x)dx A = (p f(x)dx— (p f(x)dx for % — 1, . . . , d — 1, (3.4) 

JBi JBi JB d 

and we call the number d — 1 = g of independent A- and B- cycles the genus of the 
quantum Riemann surface. 

The newly introduced A- and £>-cycles again satisfy the standard intersection for- 
mula, 

AanB p = <W, (3.5) 

and most of our construction features depend only on the homology class of the paths 
A a , B a at the asymptotic infinities, but in the intermediate considerations it is useful 
to choose a representant, the intersection point P a , 

A a nB a = {P a }. (3.6) 




Figure 4: The pattern of A- (dashed lines) and B- (dotted lines) cycles for the example 
in Fig. I3XT1 
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3.2 Filling fractions 

In random matrices, the notion of filling fractions is just the ^4-cycle integrals of the 
resolvent. Then, if the ^4-cycles are chosen to lie in the physical sheet (which is possible, 
say, in the hyperelliptic case), the discontinuity of the resolvent along the corresponding 
cuts determines the eigenvalue density and the .4-cycle integrals determine the portions 
of eigenvalues lying on the corresponding interval of eigenvalue distribution. They are 
called therefore the filling fractions. 

In the case of quantum surface, we define the "filling fractions" e a to be 



1 

2irt 



uj(x)dx = / (u)(x) — u(x))dx, a = l,...,d. (3.7) 



Aa 



Note that this definition depends on where we place the contours and (in the case 
where a sheet is bounded by more than one ^4-cycle) we also have a freedom to assign 
residues inside the sheet to different ^4-cycles, so the filling fractions are defined up to 
integers times h. 

For the difference, we have 



u(x) — u(x) 



ip a+ (x)ip a _(x) 



where w a+>a+ = 4>' a+ ip a _ — i J ' a _i J a+ is the Wronskian of the two solutions. Therefore, 
this difference decreases exponentially in sectors where the both solutions ip a+ and i/j a _ 
increase, and this is why we identify the asymptotic domains of the ^4-cycle integrals 
with "branch points." 
We have 

d 

J2^=t , (3.8) 

a=l 

which just follows from that the sum of integrals over ^4-cycles is equivalent to evaluat- 
ing the integral over Cd- This also means that we should take as independent variables 
only d — 1 = g of the variables <T a if we consider t to be an independent variable, and 
we naturally choose these g variables e a to be filling fractions corresponding to the 
cycles A a of the quantum Riemann surface. 

Remark 3.2 In the case g = —1 in j5], the only filling fraction is = to, and it is given 
by the (finite) sum of residues of the function u at the zeros sf. 

£d = to = y^Resw = h#{si}, 

I 

so to is discrete in this case. For g ^ 0, the variables € a , a = l,...,g, and to may take 
arbitrary, not necessarily integer, values. 
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3.3 First kind functions 



After defining the cycles, another important step is to define the equivalent of the 
first, second and third kind differentials. We begin with the definition of the first-kind 
differentials. 

Let hk, k = 1, . . . ,d — 1, be a basis in the complex vector space of polynomials of 
degree < d — 2. 

We introduce the functions 

^) = 7T7^T / h k (x')^ 2 a (x')dx'. (3.9) 

We use the same polynomial hk{x') for all the sheets of the Riemann surface. 

Note that because every ip a (x) is a solution to the Schrodinger equation, Vk(x) has 
double poles with no residue at the Sj (at the zeroes of ip a ), and behaves like 0(l/x 2 ) 
in the sector S a and inside all the sectors in which ip a is exponentially large (if the 
polynomial hk(x') has power less than d — 2). Therefore, the following integrals are 
well defined in the general case: 

h, a = f v k (x)dx a = 1, . . . ,g, k = 1, . . . ,d— 1. (3.10) 

If the matrix I k a with k, a = 1, . . . , d — 1 has the full rank (which we assume in what 
follows), then it is possible to choose the canonically normalized basis of h k such that 

Ik,a = $k,a- (3-11) 

The functions Vf-(x) k = 1 , . . . , g are therefore the natural analogues of canonically 
normalized holomorphic forms (1st kind differentials). 

We now extend this notion to the meromorphic (Whitham-Krichever) [15] differ- 
entials. For this, let us consider the following basis hk, k = in the space of 
polynomials of arbitrary order: the first d — 1 elements of this basis are the original 
polynomials hk each of which has degree not higher than d — 2, each polynomial hk 
with k > d — 1 has degree exactly k — 1 and must be chosen on the following grounds. 

Define the functions Vk(x) with k > d — 1 exactly as in (13.91) with hk being now a 
polynomial of arbitrary (but fixed) degree k — 1. The coefficients of hk with k > d — 1 
are unambiguously fixed by the normalization conditions: 

• (the residue condition) 

i x~ l v h {x)dx = 5 h k-d, I = 0,1,..., k > d- 1. (3.12) 
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(the normalizing condition) 

Vk(x)dx = 0, a — 1, . . . , d — 1, k — d, . . . . (3.13) 



A. 

Remark 3.3 Although the functions Vk{x) generally increase as x k ~ d as x — > oo, the in- 
tegral (|3.12p as well as the normalizing condition (|3. 13|) are well defined for any finite I and 

k. This is because the difference Vk{x) — Vk( x ) is exponentially small as x — > oo a± for any 
k, and we can integrate it along Co weighted by any polynomially growing function. The 
integral (|3.12p is therefore a natural analogue of the residue at infinity of order I + 1. 

3.4 Riemann matrix of periods 

An interesting quantity in standard algebraic geometry is the Riemann matrix of pe- 
riods provided by integrals of the holomorphic differentials over £>-cycles. So, an anal- 
ogous "quantum" Riemann period matrix r it j, i,j — l,...,g is 

r a ,i = f Vi(x)dx. (3.14) 

Note that this definition makes sense since Vi(x) (i = 1, . . . , g) behaves as 0(l/x 2 ) 
in the sectors asymptotic for the ^-cycles. And because the residues of Vi(x) vanish at 
all zeros sj , these integrals depend only on the homology class of £>-cycles. 

Like for the classical Riemann matrix of periods we have the following property: 

Theorem 3.1 The period matrix r is symmetric: ry = Tjj. 
proof: 

This result follows from Theorem 14.81 below, since: 

<£ dx (b B(x,z)dz = 2in (p dxv a (x) = 2iitTp ta 

J Bp JB a JBp 

a p p a 

and from the symmetry theorem 14.91 for the Bergman kernel, B{x,z) = B(z,x). □ 



4 Recursion kernels 

One of the key geometric objects in [11] and in [12], is the "recursion kernel" K(x,z). 
It was used in the context of matrix models [6] for constructing a solution of loop equa- 
tions. We use its analogue below for constructing the 3rd and 2nd kind differentials. 
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4.1 The recursion kernel K 



First we define the kernel 

K{%z) = T ^ r r [ X (4.1) 

SfaW Joo a X'-Z 

and for each a — 1, . . . , <?, we define 

hC a {z)=<f> K{x,z)= (K(x,z)-K( a x,z)). (4.2) 

In these expressions, we must also specify the contours of integration w.r.t. the variable 
x' from the infinities oo a± to the point x lying on the cycle A a . We assume that these 
contours go first from the corresponding infinity along the part of the adjoint cycle B a 
that lies in the sheet a± until it reaches the intersection point P a ; after this point, we 
integrate along the cycle A a towards the final point x (see Fig. 14. ip 



X £ t-^n 




Figure 5: The path of integration w.r.t. the variable x' in the expression for the 

(X. 

recursion kernel K{x,y). 



To obtain the domain of the function K(x,z), we slightly deform the contours of 
integration over edges of the *4-cycles as shown in Fig. I4.1( then, for the variable z lying 
in the domain that is "inner" w.r.t. integrations from infinities for all the functions ip 1± , 
that is, for the domain that is separated from all the infinities oo 7± by the drawn apart 
edges of the ^.-cycles, the kernel K(x, z) is well defined (and it develops logarithmic 
cuts if we push the variable z through the boundary of the sheet S a ). 

We now need to describe analytic properties of these functions. 

For a fixed x, the kernel K(x, z) is defined for z in the cross-hatched domain in 
Fig.EXU 

Taking an integration path between oo Q and x we obtain that K(x, z) is defined 
for z outside this path. Across the path }oo a , x), K(x, z) has a discontinuity w.r.t. the 
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Figure 6: The domain of variable z (crosshatched) in ( 14. ip (we slightly deform the 
^4-cycle integrals). 



argument z: 



(4.3) 



The similar statement is true for C a (z): when z cross the line of the cycle Ap, we have 

2171-^(3) f°°P± dx" 



S z Cp(z) 



h 



(4.4) 



We now define the recursion kernel K(x,z), which is the main ingredient of our 
construction. 



Oi 

Definition 4.1 The recursion kernel K(x,z) reads 



d-l 



K(x, z) = K(x, z) — Vj(x)Cj(z) 



(4.5) 



It is defined for z in the cross-hatched domain in Fig. J^.l 



Theorem 4.1 The kernel K has the following properties: 
• For a given z, 

K(x,z)~ 0(x- 2 ) (4.6) 

when x — > oo in all sectors (if the function ip a {x) increases in all the sectors 
except S a ). 
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The normalization condition reads 



(p K(x,z)dx = 0, j — 1, . . . , d — 1. 

J A, 



• K(x,z) has double poles with zero residues at the zeros of ip a . 
Theorem 4.2 We have in all sectors at infinity : 

K(x,z) ~ 0(z- d ). 

z— >oo 

More precisely we have: 

oo ,a. 

K{x,z) ~ - 2^ 



k=d-l 

with 



z k+i 



and 



K k (x) = K k (x) - ^ v j(%) f K k (x') dx'. 

7 = 1 Ja o 



proof: 

We can expand K(x, z) as 



K(x,z) ~ 



k=0 

where 



and therefore 



K k (x) = K k (x) - ^2 v a{%) f K k (x ) dx'. 

a=l ^ 

For k < d — 2, (x') k can be presented as a linear combinations of hj(x'), 

d-l 

0=1 

and from the normalization condition, we immediately obtain that 

<b K k {x) dx' = b kja , 

and therefore K k (x) = for k < d — 2, which implies that 

K(x,z) = 0{z~ d ). 

□ 
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a /3 

4.2 Third kind differential: the kernel G(x, z) 

The second important kernel to define is the equivalent of the third kind differential. 

a (3 

Definition 4.2 The kernel G(x, z) is 

G(x, % = -HiPl(z) d z = 2h^^K(x, z) - hd z K(x, z) (4-18) 

Integration by parts gives 

n r ^ 1 2 r dx' 2 (i>' a {x') %{z) 



x-z ^l(x)J 00a x'-z a \ip a (x') i)p(z) 
-hJ2 Vj (x)^(z)d z ^-. (4.19) 

j=l y^ Z > 

In what follows we often are in the situation when we take two integration contours, 
Cd x and Cd z , and must interchange the order of integration (or the order in which 
these two contours intersect the £>-cycles). From the definition of the ^4-cycles, it is 
then obvious that we must interchange the variables x and z within the same sector, 

a 

so we need permutation relations for G(x, z) with a = (3. Then, as x — > z, we find 
that G(x, z) ~ , i.e., there is a simple pole with the unit residue at z — x. Because 
the combination [ ^ a i x ,l — ^r¥\ I is regular at x' = z, interchanging the order of 

x'-Z yip a (x') tpa(z) / & & 

integration over Cd x and Cr> z then just gives the residue at z = x; no logarithmic cut 
takes place. 

Theorem 4.3 G(x,z) is an analytical function of x with a simple pole at x = z with 
residue —1 for a = (3 and with double poles at the 's (zeros ofip a (x) ) with vanishing 
residues, and possibly an essential singularity at oo. 

a /3 

G(x,z) is an analytical function of z, with a simple pole at z = x with residue +1 
for a = (3, with simple poles at z = sf \ and with a discontinuity across A 1 -cycles with 
7 = 1, . . . ,g (this discontinuity has opposite signs depending on which line of the cycle 
A~f, 7+ or7„ ; we cross; no discontinuity takes place when crossing the last cycle A n ): 

8 z G{x, P z) = ^2iitvp{x). (4.20) 



We also have 

I G(x,z)dx = 0. (4.21) 

proof: 
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All the discontinuities of K(x,z) except the one for Cj(z) (I4.4p are proportional 
to ipa(z) and vanish in (I4.18p . The discontinuity of Cj(z) gives T-27ri, and thus, the 
discontinuity of G(x,z) is ^G(a;, 2; ) = =p2i7r Vp(x). 

Since z) is regular when 2; = then it is clear that G(x, z) has simple poles 
at z = Sj, with residue —2HK(x, Sj). 

In the variable x, K(x, z) has double poles at x — without residue, and this 

a 

property holds for G(x,z) as well. 

The vanishing .4.- integral property follows immediately from ( 14. 7p . □ 



Theorem 4.4 

G(x,z) = 0(l/x 2 ) (4.22) 



when x — > 00 a /or am/ a. 

yli /arge z in the sector S 1 we have 

lim G(x, z) = G(x, oop) = r] 7i(3 t d+1 K d - X {x) (4.23) 

Z—t-OOj 

where = ±1 depending on the asymptotic behavior of ipp ~ e ±v ^ 2ti in the sheet «9 7 . 
proof: 

The large x behavior follows from theorem 14.11 The large z behavior is given by 
theorem H21 i-e. G(x,z) ~ iV'^J-fT^, 2) ~ ±t d+1 K d _ 1 (x). The sign depends on the 
behavior of the solution in this sector. □ 

a (3 

4.3 The Bergman kernel B(x,z) 

In classical algebraic geometry, the Bergman kernel is the fundamental second kind 
bi-differential, it is the derivative of the 3rd kind differential. Using the same definition 
as in [TT] , we define: 

B{x/z) = - l -d z G{xfz). (4.24) 
We call the kernel B the "quantum" Bergman kernel. 

a 

Theorem 4.5 B[x,z) is an analytical function of x. When a = (3, it has a double 
pole at x = z in the both variables x and z with no residue, and it has double poles in 
x and in z at the respective zeros Sj and Sj with vanishing residues, and possibly an 
essential singularity at 00. The discontinuity across A-cycles that was present in the 

a /3 

kernel G disappears upon differentiation, so B[x, z) is defined analytically in the whole 
complex plane. 
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proof: 

a 

These properties follow from those of G(x, z) of theorem 14.31 In particular, the 

a 

only discontinuity of G(x, z) is along the *4.-cycles, and it is independent of z, therefore 

a 

B{x,z) has no discontinuity there. □ 

4.3.1 Properties of the Bergman kernel 
Theorem 4.6 

B(x,z) = 0(l/x 2 ) (4.25) 

when x — > oo in all sectors, and 

B(x,z) =0(l/z 2 ) (4.26) 

when z — > oo in all sectors. 
proof: 

a 

Follows from the large x and z behaviors of G(x, z). □ 
Theorem 4.7 The kernel B satisfies the loop equations: 

Wa(%) J V 2{X — Z) Z J X — Z 

where P^\x, z) is a polynomial in x of degree at most d — 2, and 

H ) 0\ z ) J V 2[x — z) z J x — z 

where P% f (x, z) is a polynomial in z of degree at most d — 2. 
proof: 

a 

We begin with proving the first loop equation for B(x, z). We define: 



that is, we have 



B(x,z) = -d z ^2^-d z ) K(x,z) (1.29) 



B(x,z) = B(x,z) -^v^x) I B(x",z)dx". (4.30) 
Since (2 ^"H + d x ) Vj(x) = hj(x) is itself a polynomial of degree < d — 2, it suffices 



to prove Eq. fl4~2TD for B(x,z). 
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We have 



ip a (x) J ' ' 2 V V^O) / x- z 

i , Q V>J(*) 



(x — z) 3 ~^ z ipp(z)(x — z) ^ 
and therefore: 

^ {X) + d x ) ( B(l I) - -r±-^S + d z I MzM. I = o (4.32) 



'^ai.x) J V 2(x — z) 2 J \ X — Z 

This proves Eq. (I4.27P with 



9 p 

p(°\x,z) = -J2h j {x) <b B{x\z)dx". (4.33) 



a (3 

We now prove the second loop equation for B{x, z). We have 

<9 2 *) = - 2-^— + <9 2 <9 2 2-^— - <9 2 K(x, z), (4.34) 



where the operator J7(z) = § ^^J^j + <9 2 ) <9 2 (2^y - <9 2 ) reads 

£ (*) = ~\ dl + ^U{z)d z + ^U'(z) (4.35) 

and is therefore independent on the solution ipp(z) we started with. This is the Gelfand- 
Dikii operator [14]. We then have 



^(z) / ' i/jI(x) J OQa a V (x' - z) 4 (x' - zf x'-z / 

(4.36) 

We integrate the first term by parts three times introducing Y a (x) = ^' a (x)/ip a (x) (and 
exploiting that Y' a + Y 2 = U). The result reads 



ipp(z) J (x — z) 3 dx x — z 

This implies that 



^(-z) / ^ ' 2(x — z) 2 J dx x — z 
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*l>l(x) 7 OOQ q ' V (x' - z) 2 x'-z 

and the obtained expression is obviously a polynomial in z. The expression in the 
brackets in (I4.38j) is a skew-symmetric polynomial in x' and z of degree not higher than 
2d — 2. Moreover, all the terms with (x') k with k < d — 2 become upon integration 
linear combinations of Vj(x) and vanish identically when we apply the projection to the 
subspace of zero ^4-cycle integrals. So, the minimal power of x' that contributes to the 
answer is {x') d ~ l . But there is no term z d ~ x in the brackets because it contradicts 

the skew-symmetricity. The first nonzero term that might contribute is proportional 
to (x') d ~ l z d ~ 2 , which obviously means that the obtained polynomial p!f\x,z) has the 
maximum degree at most d — 2 in z. □ 

Theorem 4.8 We have for every a = 1, . . . , g: 

B(x,z)dx = 0, J) B{x,z)dz = (4.39) 



and 



A J A 



(t> B(x, z) dz = 2mvj(x) 



(4.40) 



proof: 

The vanishing of ^4-cycle integrals in the x variable is by construction. For the z 
variable, we have 

/ a f°°P+ a 13+ a P- 

(p B(x,z)dz = / (B(x, z — B(x, z)dz 

J Ab J oo-p 

= ~ \ G(x, oo~ + ) - G(x, oo^ + ) - G(x, oo~J + G(x, oo~_) jiAl) 

and the asymptotic conditions in all four cases for the function G in the second line 
are the same, so from theorem 14.41 we conclude that the result is zero. 
We begin with the integral over a cycle Bp. 

(p B(x,y)dy = jump of G(x, y) on Ap 

-- I G(x, ocp + ) - G(x, oop + ) - G(x, ocp_) + G(x, oop_) j 
= 2 m (l-5 M )vp(x) + 2K d _ 1 (x), (4.42) 

where we have used again the asymptotic conditions from Theorem 14.41 This formula 
implies that if we perform the integration § B for j = 1, . . . , d— 1, which is the difference 
of integrals over the cycles Bj and Bd, we obtain the formula (14.401) . □ 

The key property of the Bergman kernel is provided by the following theorem. 



22 



a P 

Theorem 4.9 B(x,z) is symmetric, 

B(x,z) = B(z,x) (143) 

proof: 

a P 

The proof uses that B(x,z) satisfies the loop equation in the two variables. We 



have: 



^p( z ) J V ^a(^) / V 2(x-z) 2 



1 2 <(^) 



[B(x, 








i/>a(x) 




x - 


■ Z 


1&(x) 




i/>a(x) 




x - 


- 2; 



^(2) y v ' - x-z j 

+ 9, ) (pf (S, *) - 9, ^ *m (AAA) 
ip a (x) J V x-z J 



We then have 



/ / / / \ N v4M_V^(f) #,(x) W 

= f 2^44 + ^ ^ ^ W ~ ( 2^ + 9x1 d x Mx) Mz) 



1pp(z) J X — Z \ 1p a (x) J X — Z 

= 2 U(x)-U(z) U'(x) + U'(z) ' 
(x — z) 2 x — z 

so that 

(x - zf (2^ + d)j P 2 (0) (x, I) + 2tf (z) + (x - 

= (x - zf (2|^| + d)j P 2 (0) (£, *) + 2C/(x) + (z- x)U'{x) 

= R(x,z). (4.46) 

Here, the first line is a polynomial in x, whereas the second line is in turn a polynomial 
in z. Therefore, R(x, z) is a polynomial in the both variables, of degree at most d in 
each variable. Moreover, we must have: 

R(x,x) = 2U(x) (4.47) 

Therefore we must have: 

R( x , z ) = L Q V'{x)V\z) - h ~ - P(x) - P(zfj +(x- z) 2 R(x, z) 

(4.48) 

where R(x, z) is a polynomial in the both variables of degree at most d — 2 in each 
variable. 
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Putting this back into (I4.46[) and using the symmetry x <->■ z, we obtain 

+ ^) (^O* J) - A (0) (Cx)) = R(x,z) - R(z,x) (4.49) 

Then, we can decompose the r.h.s into the basis hi(x)hj(z), 

d-i 

R{x,z) - R{z,x) = ^2{R id - Rj^hi{x)hj{z) (4.50) 

Applying the integral operator 

f(z) » -±-r [ Z ^(z')f(z')dz> (4.51) 

to the differential equation (14. 49 p . we obtain 

d-l 

P^(x, % - P a W(I, s) = £ (i? M - %)^(x) + A^s) (4.52) 

*.i=i 

where A\{x) is some integration constant. 

Then using the loop equations (14.71) we find by substraction that: 

( 2 W + ( B & l )- B &*)) =P?\x, P *)-P$ Si &x) (4-53) 

and applying the integral operator (14.511) w.r.t. the variable x in the sheet S a , we 
obtain 

d-i 

B(x, z) - B{z, x) = Y, ~ RjMx) Vj{z) + A(x) + A{z), (4.54) 

with A(sc) and A{z) the integration constants. 

Next, the large x and large z behavior of B imply that A{x) = A(z) = 0, and 
therefore 

B(x,z) — B(z,x) = ^~]{Rj,j — Rj,i) Vi(x)vj(z). (4.55) 
Then, using theorem 14.81 

I B{x,z)dx= I B{x,z)dz = 0V (4.56) 

J A JAj 

we obtain that 

Ri^R^ViJ, (4.57) 
which completes the proof of the symmetricity of the Bergman kernel. □ 

24 



We see therefore that our "quantum Bergman kernel" enjoys all the features of 
the standard Bergman kernel associated with a Riemann surface: it is symmetric, has 
no discontinuities, and possesses the double pole with no residue at the coinciding 
arguments (which corresponds in our case to the coinciding arguments on the same 
sheet Sk)- Using all these kernels we can then generalize the recursion of [HI El [EE] 
defining the correlation functions (see the next section). 

4.4 Meromorphic forms and the Riemann bilinear identity 

Definition 4.3 A meromorphic form TZ(x) reads 

n{x) = -±— f r(x')^ a (x')dx' (4.58) 

t*l>l(x) Joo a 

where r(x) is a rational function of x, which behaves at most like 0(x d ~ 2 ) at large x 
and whose poles r { are such that 

Res ipl(x)r(x) = 0. (4.59) 

The holomorphic forms Vj(x), the kernels G(x,z) and B(x,z) are meromorphic 
forms of x. 

A meromorphic form 1Z(x) defined by (I4.58P has poles at x = r iy the poles of r(x), 
with degree one less than that of r(x), it has double poles at the with vanishing 
residues, and it behaves like 0(x~ 2 ) in all sectors (having also an accumulation of 
poles along the half-lines Li of accumulations of zeroes of ip a ). Note that the integrals 
§ A 1Z(x)dx are well defined. 

We then have the following theorem. 

Theorem 4.10 Riemann bilinear identity 

Consider a meromorphic form TZ(z). For z in the sector S a (outside the 



crosshatched domain in Fig. |^.1[), we obtain the representation formula 



U(z) = -J2J2 ResG ^>^ n & dy -Yl Yl ^sG(z,y)TZ(y)dy 



p ri&Sp p m 



K 



+ y>i(l) / n(y)dy. (4.60) 



proof: 

We begin with the integral over Co (Fig- 13. 1 p of G(z, y)1Z{y)dy: 



/ G(z,y)K(y)dy = J2 G(z,y)K(y)dy. (4.61) 
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a 

This integral is identically zero because of asymptotic conditions G(z,y) — > 1/y as 

y — > oop and lZ(y) ~ 1/y 2 as ?/ — >• oo^ and by the fact that no accumulation of zeros 

occurs for the function 1Z(y) on the boundaries between the sectors Sp and Sp±i. We 
then push the integration contours through the complex plane towards the .A-cycles 
as in Fig. 13. 1. H the residues at the points 7™j and s£ give the first line of (I4.60p . the 
residue at the point x — y in the sector S a gives the left-hand side because we have 
from Theorem 14.31 that G(z,y) — l/(z — y) + reg.; it remains to consider the integrals 
along the ^4-cycles. For this, note that our integrations over y are outer w.r.t. the 
integrations over the variable x along ^4-cycles in the formula (14. 2p for the factors Cj, 
and when we push the integration over y through that for x, we have the discontinuity 

of G(z, y), which is equal to 2iTiVj(z); no such discontinuity occurs for the cycle Ad- 
All these discontinuities are independent of y, so the contour integral of the product 
factorizes for each cycle Aj and gives the second line of f)4.60p . 

To evaluate the remaining integrals inside the crosshatched domain in Fig. l4.1l recall 
that G(x,y) = ipp(y) d y (^K(x,y)/^(y)J, so integrating by parts we obtain 

( a p + £ + a p p \ f°° 



G(x. z YRi : )-G(x. : )K( :))«/: = - / {K(.v. z)r(z)-K(x, z)r(z))dz = 0, 

(4.62) 

and these contributions vanish for all the cycles Ap. □ 



5 Correlation functions. Diagrammatic represen- 
tation 

In this section, we define the sectorwise defined versions of the quantum correlation 
functions from paper I (deformations of "classical" correlation functions introduced in 
[9J El [12]). Our definitions are inspired from (non-Hermitian) eigenvalue models (see 
section |8]), but they are valid as well in a general setting of an arbitrary Schrodinger 
equation. 

5.1 Definition and properties of correlation functions 

Definition 5.1 We define the functions Wn (xi,..., x n ) called the n-point correla- 
tion function of "genus" h by the recursion 

W{ ' (x) = u(x), W£'(xi, x 2 ) = B^x-l, x 2 ) (5.1) 
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W&CSo, J)=<f dx K( a x°o, x) (w^\x, x, J)+J2 E W$ +1 (x, I)W^ +1 {x, J/ 1)) 

JC-Dx r =0 IcJ 

(5.2) 

where J and I are the collective notation for the variables (J = {x%, . . . , x n }), the 
symbol means that we exclude the terms (r = 0,1 = 0), (r = 0,1 = {x^}), 

(r = h, I = J/{xi}), and (r = h, I = J), the integration over the contour Cd x is 
defined in ( Iff. and where 

W y n \x 1 ,...,x n ) = WW(x 1 ,...,x n ) - -j — ■ . 5.3 

2(xi - x 2 y 

Here the point x is outside the integration contour Cd x for the variable x and all the 
Xi are outside the A-cycles of the projection integrals. 

To shorten equations we introduce the notation 

Ui h \x, J) = Wt?& 2, J) + EE<i^ nWt^ix, J 1 1)- (5-4) 

r=0 IC.J 

The main property of 15.11 is that these quantities solve the loop equations in the 
l/iV 2 -expansion. We also prove the following properties: 

Theorem 5.1 Each Wn h \%i, . . . , x n ) with 2 — 2h — n < is an analytical functions 
of all its arguments with poles only when x^ — > s^ l \ It vanishes at least as O (1/xf) 
when Xi — > oo ai and has no discontinuities across A-cycles. 

Corollary 5.1 We have that 

/ W^ x {x, J)dx = t 5 nfi 5 K0 . 
J c Dx 

proof: 

We proceed by recursion on 2h + n. The analyticity is obvious; the theorem is true 
for W 2 \ Assume it is true up to 2g + n, we shall prove it for W^I x {xq,X\, . . . ,x n ). 
To prove the asymptotic behavior, note that the definition 14.11 implies, first, that the 
term J Cd K(x,y)Un h \y, J)dy is of order of l/(ijj 2 (x)) ip^(x')dx' /(x') 2 ~ x~ d ~ 2 , so 

the leading contribution comes from the terms proportional to Vj(x) ~ x~ 2 , which 
completes the proof. □ 

We also have the following simple lemma that follows from the corollary 15.11 and 

ex. 

from the normalization conditions for the kernel K(x,y). 
Lemma 5.1 For all (n,h) ^ (0,0) we have 

I W^ 1 (x,J)dx = 0. (5.5) 
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Now comes first of the main theorems. 

Theorem 5.2 For 2 — 2h — n<0, the satisfy the loop equation. This means that 
the quantity 



( W ( n h \x, J/{x 3 })5 a , a] - wf&v J/{x s }) Y 

zs a polynomial in the variable x, which is of degree at most d — 2 and is independent 
on the choice of the sector S a . 

proof: 

in appendix [S] □ 

Theorem 5.3 Each wi h) is a symmetric function of all its arguments. 
proof: 

in appendix [Cj with the special case of in appendix [B] □ 

Theorem 5.4 For 2 - 2h - n < 0, Wk \ Xi, . . . , x n ) is homogeneous of degree 2 — 
2h - n: 

j=i J i=i / 
= (2-2/ i -n)^(x 1 1 ,..., Q x n ). (5.7j 

proof: 

Under a change — >■ At/-, h — >■ A/?., e.j — >■ Ae.j, and t — ► Aio the Schrodinger equation 
remains unchanged, and thus ip is unchanged. The kernel K is changed to K/X and 
nothing else is changed. By recursion, Wri is then changed by \ 2 - 2h ~ n . □ 

5.2 Diagrammatic representation 

The diagrammatic representation for the correlation functions structurally coincides 
with the one for the correlation functions in one- and two-matrix models [9j HI |6]. We 
introduce the three kinds of propagators 

x ^ y x y x y 

K(x,y) G(x,y) B(x,y) 
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and assume the partial ordering from "infinity" to "^4-cycles" to be from left to the 
right in graphical expressions. We represent the terms Wn\j) via the graphs with 
three-valent vertices; we assign its own variable £ to every inner vertex and assume the 
integration w.r.t. this variable along the contour Cd, the order of integration depends 
on which vertex is closer to the ".A-cycles" : we begin with integrating at the innermost 
vertex. We also have n outer legs (one-valent vertices) corresponding to the points x i: 
i = 1, ... ,n. They are assumed to lie outside all the inner integrations. For example, 
the term W^°\xi, X2, £3) then has the form 





whereas the recurrent relation (15. 2p takes the form 




We now formulate the diagrammatic technique for constructing the functions 
wi h \j) for n > and 2g — 2 + n > 0. Formally it is the same as the one in [9j [6]. 

We comprise all the diagrams with the corresponding automorphism multipliers 
such that 

• for W$p(J) a dia gram contains exactly n external legs and h loops; 

• we segregate one variable, say, x±, and take all the maximum connected rooted 
subtrees starting at the vertex x\ and not going to any other external leg; 

• we associate the directed propagators K(x, y) with all the edges of the rooted 
subtree; the direction is always from the root to branches; 

• all other propagators that comprise exactly h inner propagators and n — 1 re- 
maining external legs are B(x,y) if the vertices x and y are distinct and B(x,x) 
for the loop composed from the single propagator; 
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• each rooted subtree establishes the partial ordering on the set of three-valent 
vertices of the diagram; we allow the inner propagators B(x,y) to connect only 
the comparable vertices (a vertex is comparable to itself). 



6 Deformations 

In this section, we consider the variations of correlation functions Wn 9 ^ under infinites- 
imal variations of the Schrodinger potential U(x) or h. Infinitesimal variations of 
the resolvent u(x) can be decomposed on the basis of "meromorphic forms" Vk(x), 
k — 1, . . .; we set these forms to be dual to special cycles with the duality kernel being 
the Bergman kernel. It turns out that the classical h = formulas retain their form 
for H ^ 0. 

6.1 Variation of the resolvent 

We consider an infinitesimal polynomial variation: 

u->u + 8u, h^h + 5h. 

Since U = V' 2 /A - hV" /2 - P, we have 

Su=—5V'--5V"-—5V"-6P. (6.1) 
2 2 2 K ' 

We can consider variations of V'(z) w.r.t. the higher times tk, k — 1, . . ., as well: 

5V'{x) = ^5t k x k -\ (6.2) 
k=\ 

Then, for k < d + 1 , 5P is of degree at most d—1 whereas for k > d + 1 the polynomial 
5P has the degree at most k — 2. 

Computing 5 (tp' a (x) /ip a (x)) , we obtain 

5 mx)/U*)) = £ (W) - 2 f £/(*')) dx>, (6.3) 

and for u{x) = V f (x)/2 + hip' a (x)/ip a (x) we have 

6.2 Variations w.r.t. "flat" coordinates 

We choose a system of "flat" coordinates on our genus- (d — 1) manifold: 

ei, • • • , €d-i, h, h, ■ ■ ■ ■ (6.5) 
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6.2.1 Variations w.r.t. the filling fractions 

For the filling fraction 5e a we have 5V = and thus 

5U(x) = -5P(x) (6.6) 
where deg8P < d — 2, so we can decompose it in the basis of h a : 

5P(x) = J2 c »'h a '. (6.7) 



So, from (16.41) . we have 

8u(x) — — c a ' Vqt(x), (6.8) 

a' 

and because 2me a i — §^ f cj, we obtain 

2«7r<5 a , Q ' = ® 5w = — ® c a » f a » = —c a i (6.9) 



Therefore 5U(x)/5 Ea = 2mh a (x) and 



5 ta io(x) = 2mv a (x) = (p B(x,z)dz. (6.10) 



The flat coordinate e a is dual to the holomorphic form v a , which is itself dual to 
the cycle B a : 



e a = — j> w, S €a u = 2mv a = I B. (6.11) 



6.2.2 Variation w.r.t. to 

We have 

SU(x) = -5P(x) = -t d+1 x*' 1 + Q(x), (6.12) 
where degQ < d — 2. Using Eq. ( 16. 4 p we obtain 



8u{x) = -±- f (-t d+1 x' d - 1 + Q{x')) ^ a {x') dx', (6.13) 

€W Joo a 

and the polynomial Q must be chosen such that § A 5u> = We therefore have 
5uj(x) = -t d+1 K d _x(x) 

d-l „ 

= -t d+1 (K^tix) - y2 v p(?) f K d _ l (x')dx 



/3=1 



^(a), (6.14) 
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where 

Kk(x) = 3^ f X x ' k dx, (6.15) 



and Kk(x) is the kth term in the large z expansion of K(x, z) = — ^2T=o ^ife^ com- 
puted in theorem 14.21 From theorem 14.41 we have G(x, oo a ) = r] a td+i Kd-i(x). This 
shows that 

a 1 a a r°°d a 

S to u(x) = - (G(x, oo d - + ) - G(x, oo 3 J) = / + B(x,z)dz. (6.16) 

The integral is taken here over the last, dth I?-cycle. 

The flat coordinate to is then dual to the 3rd kind meromorphic form — 2G(x, oo.), 
which is itself dual to the cycle [ooj_, ooj ]: 



t = f oo(z)dz, 8 to uj(x) = / B(x,z)dz. (6.17) 

JC D J oo j_ 

6.2.3 Variations w.r.t. t k , k = 1 . . .. The two-point correlation function. 

Because 



t k = <b htp'(z)/tp(z)z- k dz, k = 0,l,..., 
Jc D 

the conditions dtk/dt r = 5 fc r and dt k jde^ = imply 

and from the general form (I6.3P of variation, we conclude that (cf. ( 13 . 12j) and (I3.13P ) 

i.f#U^). (6.19) 



8t r \ 1pa(x) 

In turn, we have the following lemma. 
Lemma 6.1 

a If a Z T 

v d+r (x) = — (h B(x,z)—dz, r = l,.... (6.20) 
2m Jc D >x r 

proof: 

That the expression in (I6.20p has the desired structure follows from the explicit 
form of the kernel B\ we need only to verify the normalization conditions. Vanishing 
of „4-cycle integrals is obvious; we must verify 

$d,i — 77~ (f x~ l Vd+ r {x)dx = . , . (f dz (f dx — B(z,x)x~ l . (6-21) 



JCn ""' " ( 2m r JC D >x JC D 
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Interchanging the order of integration contours and using that x~ l B(z,x) ~ x~ l ~ 2 as 
x — > oo, we observe that the only nonzero contribution comes from the double pole at 
x — z, which gives 



dxx r 1 1 = Sri- □ 



2wr f Co dXX (dz r ) Z=;E 2m f Co 
□ 

We now define the operator of the loop insertion 

_d_ = f. 

applying which to hip' /ip, we obtain 

and since j> Co B(x, z)dz = 0, we add the term with r = into the sum obtaining 

a 1 

B(x, z) dz 

y>C D y - z 

in the right-hand side. Note that the point y lies between some infinity, say, oo^ and 
the integration contour Cd- Pulling the contour of integration through the point y 

a /3 

to infinity we obtain zero due to the asymptotic conditions for B(x,z), so the only 
nonvanishing contribution comes from the residue at y — z, which eventually gives 

1 h-tfj-!- = — B(x,y). (6.23) 



dV(y) V Mx)J 2 
Correspondingly, since V'(x) = ^1^2 , we obtain for the two-point correlation 
function W^\x^y): 

6.3 Variation of higher correlation functions 

Note that for all the above variations w.r.t. the flat coordinates, we have a cycle 5lu* 
and a (sector-independent) function A* SuJ such that 

6u(x)= f B(x,z)A* Sul (z)dz. (6.25) 

Jslu* 

The theorem below allows computing infinitesimal variations of any under a 
variation of the Schrodinger equation. 
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Theorem 6.1 Under an infinitesimal deformation U — >• U + 8U, we have: 

5W^(x\, . . . , a x n ) = [ Wi%(x\, £», aO A*(aO dr* fS.itfJ 
where (Su*, Ag u ) zs i/ie ditaZ cycZe to the deformation of the resolvent oj — > u + Su. 
proof: 

We prove this theorem by induction. We begin with the loop equation for 

Wi h) (x,j): 

(^Wl + hd^j W^ix- J) + Ui h \x,x; J) = PW(X, J). (6.27) 

Taking a variation 5 w.r.t. any of the flat coordinate, we have 

^tM^ + m) 5W^\x,J)+(25h^) Wi h \lj) + 5U( h \x,x;J) 
= SP( h \x,J), (6.28) 

where 5P^ (x, J) is a polynomial in x of degree at most d — 2. Here both 8U^ {%, x; J) 
and 26h%$ = J^. 2B (x,x')A*(x')dx' can be expressed by the induction assumption 
in the dual-cycle-integration form; meanwhile 

5Ui h \x,x;J) = ! U^ix, x; J, x')A*(x')dx' - [ 2B(x, x')A*(x')dx' ■ W^ h \x, J), 



8u>* J Slu* 

(6.29) 

because no term containing the two-point correlation function W^\x,x') enters 
5U^ h \x } x;J). 

Using the loop equation of form (I6.27P relating J,x') and U^^x, x; J,x'), 

we observe that the second term in the r.h.s. of (I6.29p cancels the contribution of 
28h t a ¥\ , and we obtain that 

2h^f\ + Hd x ) ( [ (x, J, x')A* (x')dx' - SWP (x, J) 



= SP( h \x,J)- P}$i(x, J,x')A*(x')dx' 

d-1 

= J2 a i( j ) h i( x ) ( 6 - 3 °) 

1=1 

with the r.h.s. being a polynomial in x of degree not higher than d — 2 expressed in 
the basis of the polynomials hi{x). Using I6~4"t we have 

» d-l 

/ W^%(x, J,x')A*(x')dx' - 5Wi h \x, J) = Y^ a i( J ) v i( x ) ( 6 - 31 ) 
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but since both wi h \x, J) and W^+i(x, J,x') have vanishing ^4-cycle integrals, we have 
that ctj = 0, i.e. 

5Wi h) (x,J)= ! W^ 1 {x,J,x')h*(x')dx l (6.32) 

Jul* 

□ 



Corollary 6.1 



dwi h \.J) 



f ^%(J, x')dx' for any n>0,h> 0. 

JBr, 



7 Classical and quantum geometry: summary 



In the table below we summarize the comparison between items in classical algebraic 
geometry and their quantum counterparts. 





classical h = 


quantum 


plane curve: 


/•;(,-.//) /•:,,;.'••//■' 

E(x,y) = 


e(x, y) = E i,3 xi y j > b> ^1 = ^ 

E(x,hd x )ij = 


hyperelliptical 
plane curve: 


y 2 = U(x) 
degU = 2d 


- [/ (x) , [y, x] = h, 


Potential: 


V'(x)=2(y/U(x)) + \ 


resolvent: 


u(x) = V'(x)/2 + y. 


uj(x) = K'(x)/2 + | 


physical sheet (s): 


y ~oo -^V(x), UJ ~ t /x 


v 

sectors where ip ~ e~2i 


branchpoints: 


simple zeroes of U (x) 
U( ai ) = 0, U'( ai ) ± 
i = l,...,2d + 2 


half-lines of accumulations 

of zeroes of ^ 

L t , i = l,...,2d + 2 


double points: 


double zeroes of U(x) 
U{ai) = 0, U'(ai) = 


half-lines without accumulations 
of zeroes of ^ (?) 


genus g — — 1 


degenerate surface 


^ e v/2fi =p lynomial 
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classical h — quantum 


A. a -cycles 


surround pairs of surround pairs of half-lines 


01 = 1 a 


branchpoints 1 of accumulating zeroes 


extra .Ad-cycle 


surrounds last pair of 1 surrounds last pair of half-lines 


a — d 


branchpoints of accumulating zeroes 


£>-cycles 


At n Bj = 6ij 


Holomorphic 


vAx) = ~ h ^ x \ ' Vi(x) = , s r ipl(x')hi(x')dx' 


forms, 1st kind 


hi =polynomials, deg hi < d — 2 


differentials 


normalized: £ . vAx) dx — 5 a a — 1, . . . , d — 1 


Period matrix 


T~i,j — §q. Vi , i, j = 1, • • • , Q , Tjj = Tjj 


Filling fractions 


2m e a = § Aa u{x)dx , a = !,...,$, e d = t - E«=i e " 


3rd kind form 


G(x,z) ~ x ^ z \/{z-x) 
G(x,z) = (2u(z) - V\z) - hd z )K(x,z) 


Recursion kernel 


K(x,z) 

K(x, z) = K(x, z) - J2i Vi(x) Ci(z) 
Ci(z) = § A _ K(x',z)dx' 

K(X, Z) — , K (X, Z) — i ,$ , \ i 

V ' / z-x 2^/U{x) 1 V ' Wiw ™o x ~z 


Bergman kernel 


B(x, z) = — i a z G(x, z) 


r\ i i * i i • rr j * i 

2nd kind differential 


T>/ a a \ 1 In/ \2 

B{x, z) ~ \/2{x - zy 


Symmetry: 


D\X^Z) = jDyZ^X) 




§ A B(x, z)dx = 




§ B _ B(x, z)dx = 2iirvi(z) 


Meromorphic 


n ^ dx ~ 0lm I n & ~ r ^ ^* V) dx' 


forms 


r(x) =rational with poles r iy r(x) = 0(x d ~ 2 ) 
| Res ri r(x')tp 2 (x') = 


Higher 


Wi%(x, J) = Y. i2 h § Ci K & *) ^ (W^iz, z, J) 


correlators 


+ £',+,' h twt> i W H\ T \ (z, i)w\ s 2 m (z, ry) 

Z_-/S+S — fi, IWl — J l + |l|V ' / l+|i'|V ' / J 

where Ci surrounds the branchpoint lu^ 7^— <£, = <£, 


Symmetry 


h4 9, (xi,x 2 , ...,x n ) = W^\x u{l) ,x u{2 ), ■ ■ -,x a(n) ) , a e S n 


Variations and 


U(x) ->• + 


dual cycle 


5C/*: 6u(x) = f su ,B(x,x')A su {x') dx' 




hM*) = in B ^' x 'y £ vdx' 


variation 5t 


S to u{x) = J°° s + B(x,x')dx' 

d_ 


variation 5t{ 


5 ei cu(x) = § B B(x, x') dx' 


Variations of 
higher correlators 


SW^ix-i, ...,x n ) = j 5u , Wj+^xi, . . . , x n , x') Asu(x') dx' 
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8 Application: Matrix models 

The main reason for introducing is that they satisfy the loop equations for the 
random /3-eigenvalue ensembles. We can therefore identify them with the correlation 
functions (resolvents) of these ensembles. 

Consider a (possibly formal) matrix integral: 

Z= [ rfMe-^ trV(M) (8.1) 



En, is 

where V(x) is some polynomial, and where E^^ = is the set of hermitian matrices 
of size N, £>jv,i/2 is the set of real symmetric matrices of size N and £7jv,2 is the set of 
quaternion self dual matrices of size N (see [TT]). 

Alternatively, we can integrate over the angular part and get an integral over eigen- 
values only [T7] : 

r N 
Z = d\!...d\ N |A(A)| 2/3 n e ~^ y(A,) ' ( 8 ' 2 ) 

i=i 

where A (A) = rii<j(^j ~ ^t) * s the Vandermonde determinant. 

We can then generalize the matrix model to arbitrary values of (3 taking the integral 
(18. 2 p as a definition of the /3-model integral. For this, we take 

n=%(vp-Vi- (8-3) 



Notice that h = correspond to the Hermitian case (3 = 1, and K — > —h corresponds 
to p -> 1/(3. 



8.1 Correlation functions and loop equations 

We define the correlation functions (the resolvents) 

W k (x 1 ,...,x k )=(3 k / 2 ( V _J_...^L_\ (s.l) 

\ . ^ x x - X h x k - X ik / 

and 



W = 7 = hgZ. (8.5) 

When considering variations in the potential V(x), we again assume these resolvents 
to satisfy the asymptotic conditions sectorwise, which means that they are defined also 
sectorwise. 

And we assume (this is automatically true if we are considering formal matrix 
integrals), that there is a large N expansion of the type (where we assume h = 0(1)): 

oo 

W k (x\, . . . , x\) = J2(N/to) 2 - 2h - k wi h \ a x\, . . . , a x\) (8.6) 

h=0 
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W = F = ^(N/tof-^W™ = ^(iVAo) 2 " 2 ^. (8.7) 

h=0 h=0 

The loop equations are obtained by integration by parts, for example: 



gives: 



(8.9) 

We define the polynomial 



v / i 2 o V i i iV^m) 

^\(x-A,) 2 + '^-.r A, A, A, t x->. 
^-r 4 \{x — Xi) 2 x — Xi x — Xj t x — X,i 

Y I h/3 y 1 1 AV5m) 

■^-^ \ (x — Aj) 2 ^-r 1 x — X-i x — Xj to x — Xi 

09 - + /?(^! 2 (x) + iwaCn, *)) 



P l(l) = ^ E (m^(M) = ( K^ 1 ) + . (8^0) 
We thus have the loop equation of [3] 

W?(x) + h— W[{x) + W 2 (x, x) = — {V'{x)W l {x) - Pt{x)) (8.11) 
Using the expansion (18 .6p we come to the Ricatti equation 

W[°\xf + hd x Wi°\x) = V\x)W[°\x) - pf\x) (8.12) 

satisfied by u(x): 

Wl°\x) =u(x). (8.13) 

The correlation functions of /3-eigenvalue models obey therefore the topological recur- 
sion of definition 15.11 
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8.2 Variation w.r.t. h 



In this subsection, we use the analogy with the /3-eigenvalue ensemble to hint the 
possible form of the last remaining building block of our construction, which is the 
variation w.r.t. h, the exponent of the Vandermonde determinant in (|8.2|) . Up to 
irrelevant multipliers, we can consider (3(d/d(3) instead of h(d/dh), for which we have 

/3-^logZ ~ ^ J d\ 1 ...d\ N A(A) 2 ^log|A(A)in e -^ y(Al) , (8.14) 

so the logarithm of the Vandermonde determinant appears. 

Ct 

It seems impossible to construct such a term from Wi(x), but we can use the two- 
ex 7 

point correlation function Wzix^y) instead. Adopting a /3-model inspired definition of 

a P 

W2(x,U) as a two-resolvent correlation function (not necessarily connected), 

1 r N 1 N 1 m 

J 1=1 * j=\ y 3 

we then introduce the regularization (both IR and UV, if speaking in the physical 
terms). At this point, we also split all the eigenvalues A, into clusters each of which 
corresponds to some sector 5* 7 ; for each term l/(y — Aj) we then integrate w.r.t. y 
from A 7 to x + <5 7 along the straight lines all of which are parallel; the regularization 
parameters depend only on the sector number 7, and the limit of removed regularization 
corresponds to A 7 — > oo 7 and 5 7 — > 0. We then obtain 



2/3 f X+S -> a 7 



J i=l % 7 ?-y=1 A ^ j7 i=l 

TV 

= / dx 1 ...dx N A{x) 2 ^Y,—r 
J i=i x ~ 

x E E [ lo § l x + - A *rl " lo S 1^1 + °(VA 7 )] II e~^ VM . (8.15) 
7 i7 =1 i=1 

We now want to perform the integration over the variable x to obtain the expression 
(I8.14p . Obviously, we need to choose the integration contour in a rather specific way: 
we want it to encircle all the poles in Aj in the variable x leaving outside all the 
logarithmic cuts from oo 7 to Aj + 5 7 in the corresponding sector (see Fig. 18. 2p . Given 
such a contour, we can then perform the integration w.r.t. x by residues at the points Aj 
(recall that in the eigenvalue model pattern we do not have boundaries between sectors 



e s o 



7 j~t=l ' 

N 

1 



N 
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inside the complex plane as yet; they appear because of the collective effect of taking 
into account the A poles and due to sectorwise regularization chosen). Evaluating the 
integral over x in (I8.15P by residues at Aj, we obtain 

N e 7 

dAi . . . dX N A(A) 2/3 J2 J2 log I Xi + 5 -y ~ Xi 



20 
Z 



\?7 I 



»=1 7 3~t =1 
N 



-iV]Te 7 log |A 7 |+0(1/A 7 ) \l[e 



«0 



V(Xi) 



i=l 



^ | dX 1 ...d\ N A(\r 



log JJ(Ai - Xj + 5 n 



#3 



-2PN ^ log |A 7 | + 2/3 ^ e 7 log |5 7 |, (8.16) 

7 7 

where the first term in the r.h.s. gives the desired integral (I8.14p as 5 1 — > 0, and the 
last two terms are divergent in the limit of regularization removed. These two terms 
depend however only on the occupation numbers therefore contributing to potential- 
independent part of J-q only and we can remove them by the proper normalization. 




Figure 7: The origin of the integration contour Cd in the matrix-model concept. The 
inner dots are Aj and the outer dots are Xi + <5 7 (7 = 0, 2, 4, 6); thin arrowed lines are 
the logarithmic cuts. 



9 The free energy 

We use the variations and theorem 15.41 to define the T^. 
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9.1 The operator H 

Theorem 15.41 gives: 



d+l 



(2-2h-n- hd h ) = [t d to + E ** \ + E ^ ^ (9-1) 



fe=i 



In section El we expressed the derivatives of Wn as integrals of W!£l 1 up to the action 



i=l 



Of^, 



(2 _ 2A - n - ftd ft ) WW = ff. = H.-^j W$* 



d 



dV 



(9.2) 



where H is the linear operator acting as follows: 



H.f(x)=t I f + J2 [ '-^f • » # 



(9.3) 



Bi 



We define = for n = and h > 2 as 



Definition 9.1 TTie /ree energy J^h for h > 2 is the functions for which 

(2-2h-h d h ) T h = H. W[ h) (9.4) 

9.2 The derivative H-j^ 

The matrix-model considerations in preceding section imply that constructing the 
derivative in H of the correlation function Wn (J) would involve resolvents of order 
n + 2. That is, 



„_n rr- r •/ OO 



r=0 JCJ 



(9.5) 



where £ must be taken to be an "innermost" variable in the sense that taking into 
account that j B(£',y) = G(£,y), we replace all the appearances 
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and 





with no additional factors. 

Note that the sum in (19 .5p ranges all cases, not necessarily stable ones, so we begin 
with studying nonstable contributions to stable cases (2h — 2 + n > 0). Note that all 
these contributions then come from the second term in ( 19. 5ft . 



9.2.1 Case r = 0, I 



and r = h, I = J 



We consider the situation when n > 1. We can then fix x\ to be the root of all the 
subtrees composed from the .fT-propagators and £ can then be the variable of any of 
external S-legs. Then, the contribution in W^\j) comprises all the insertions 



+ 



(o) 



and 



(o) £ £ 
+ ~^rf f 



(0) 



Let us consider the first case; the second one can be treated analogously. We push the 
integration over £ through the one over r\ in the second diagram and as the result we 
obtain 



7] 7] 1] 

Here the sum of the first two terms contains the integral of the total derivative of 
the function fW^ffldgGfat), and since G(r/,£) ~ O^ 1 ) and f W{ 0) ~ 
J^tod^'/^' ~ to log |C| 3 this contribution vanishes. Only the third contribution coming 
from the residue at £ = rj survives, and this contribution is nothing but minus the 
action of the H operator on the external leg .8(77, £), so 



" - ** w =0, 



V 



So, the total contribution of the two cases r = 0, I = and r = h, I = J exactly 
cancels the action of the H operator. 



42 



9.2.2 Case r = 0, I = {x x } 
We begin with the identity 

G(x,t)K(£,y) = -K@,y). (9.6) 

>c D /. >y 

(Here and hereafter inequalities of the type x > Cd 6 > y indicates the mutual positions 

of points and integration contours.) Indeed, representing G(x,£) = ^(O^g and 
integrating by parts, we obtain 



T+ 



x>C Dj _>y 



Z-y 



+ X>i(£)GKv) 



the substitution apparently gives zero, and in the second term only the residue at 
£ = y contributes thus producing (19.61) . An obvious corollary is the second convolution 
formula 

r ™ 8 ~ b 

(9.7) 



G(x J 0B(d) = -B(x,y). 



x>C D( ,>y 



In the case r = 0, I = {xi}, we have the diagram 



Wk h \^j/{x x }) 



wi h \j) 



9.2.3 Case r = h, I = J/{x n } 
Here, we need another identity 

/ G(x,t)B{t,y) = 0. (9.8) 

Jx,y>C D ^ 

to obtain it, we represent the functions G and B through the kernel K, that is, we 
have 

/ G(x,£)B(d) = 
Here the substitution gives zero, and the third-order differential operator acting on the 

8 

kernel K(y,£) is again the Gelfand-Dikii operator (I4.35P that does not depend on the 



43 



sector 7. The integrand is also obviously regular at all zeros of ^-functions, so the total 
integration over ^4-cycles just gives zero. 

Therefore, the contribution of the case r = h, I = J/{x n } is zero, and the total 
contribution of all the unstable cases together with the action of the if-operator just 
gives minus one times the original contribution Wn(J). 

9.3 Examples of application of hJ^ 

9.3.1 Reconstructing wi° ] (J) 

We now apply formula ( 19. 5 j) to reconstruct the correlation function Wn\j). In the zero 
genus case, we must take into account only the nonconnected contributions (the second 
term in (I9.5P ) into account; we choose the root of the first term, W$ +1 (£' , I)d£', 
to be x\, the point £ is then the end of some other (nonrooted) leg G(r),£) of the 
first diagram; for the second diagram we choose the root to be at the end £ of leg 
with the corresponding propagator K(£,p). As the result of the integration over £, 
we obtain using (I9.6P that these two diagrams are sewed along the propagator K (77, p) 
thus producing the connected diagram with the maximum subtree of propagators K 
rooted at the external point x\. We may now ask the question how many times the 
given diagram can be obtained as a composition of two diagrams in the formula (19. 5p ? 
We obtain this diagram by first breaking it into two parts by cutting some of internal 
arrowed lines (including also the external line K(xi,k) if we take into account the 
nonstable contributions) and then sewing again along the same line; obviously, we 
obtain this diagram as many times as the total number of arrowed lines (with the 
minus sign from (19. 6p ). which is 2 — n for Wn\j). So, we see that adopting the 
definition ( 19. 5 p for the action of the operator ftjjr, we obtain 

{ h ~k + s -w) W - ){J) = {2 ~ n ^ ] ^ ^ 

which is a particular case of formula (19.21) . 

9.3.2 Acting on W 2 \x u x 2 ) 

Here, we consider the action on a nonstable correlation function W 2 {xii ^2) — 

olx cx.2 S 

B{Xi, x 2 ) — ( g "Iw)2 ■ Excluding the terms that compensate the action of H, we have 
that the action of h-j^ gives 

dtG@i,0(B(Z,7 2 )- „ 1 „ )- [ d£—!—B($ 2 ,t) + x 1 <r>x 2 



(£ — X 2 ) 2 J Jxi,X2>C D * X 1~Z, 

2 °i \ _ 
: 2; X l> — 

which again is in accordance with formula (19. 2p . 



X1>Cd£>%2 " -^2) J J xi,X2>Co f 

-2BCx\, x 2 ) + B(x\, x 2 ) + B( a x\, x\) = 0, (9.10) 
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9.3.3 Acting on T\ 



Here we demonstrate the first case of "reconstructing" the free energy term (although 
this is not true, to obtain the genuine term T\ we need other methods, which are still 
missing; we must however demonstrate that applying formula (I9.2p we get zero perhaps 
up to some irrelevant regularizing factors). In this case, no nonstable terms contribute; 
the only contribution comes from the first term in (19. 5p . which gives 



h dh +H -W 



di G{t, i) 



?+5 >' d ^')/MO-i 



i-i 



di'h 3 {i')M0) 



MO 

mo 



Integrating by parts in the second term, we obtain (up to terms of order 0(5 a )) 
Ic D d£ hjiQCjiQi an d the integrand is sector-independent and nonsingular at zeros of 
if) a thus giving zero upon integration. In the first term, integrating by parts in the vari- 
able £' the term with l/(£ — £') 2 ) and taking into account that lim^/ 



1 



we obtain 



di 



Cii, 



di 



di 



di 



di 



di 



di 



MO 
MO 



di' 



+ 



MO 

MO 
. MO 

MO 
MO 

MO | 
MO 

' MO , 



x+S a 
x+<5 Q 



2 <(Wa(0 



c-e mo 

2 <(w«(^ 
c-e €(0 
M0)MZ 



+ 



g(g) 



MO 

i 



e-e mo 



di' + 



d- 



9Me) de+ MQ d 



Ma 

MO 

i 



MO e-e 



^(0 e-tMO 



+ 0(5 a 



x+S a 



and the result is a constant, which is divergent in the limit of regularization removed 
but is otherwise independent on all the variables (the same phenomenon occurs when 
calculating the corresponding action of the H operator on T\ in the standard matrix 
models, see jU [3]). 
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9.3.4 Acting on W{ 



Mi 



In this case, we have two possible contributions: the one from nonstable graphs gives 
Wi (x) with the (desired) factor —1 whereas the second one would come from the first 
term in (19.51) originated from term, that is, 



dr)K(x,r))G( V ,OB( V ,0, 



(9.11) 



where the contour of integration over rj goes between the points £ and £. We set 
integrals of this type to be zero, which provides the last required prescription for the 
diagrammatic technique describing the free energy terms Th- 



9.4 The term T h 

For the stable cases (h ^ 0, 1), we can now formulate the diagrammatic technique for 
the term Th- We need the diagrams describing the stable terms Wi T \£) and W^ h ~ r \^) 



with 1 < r < h — 1 and W. 



(h-i) 



(CO- 



(2h - 2)F h 




- wt 1] 77 

p 1 


h-1 I 








-E ^ 

r=l 






rj 1 



Here in the first term the sum ranges all the diagrams contributing to that 
have distinct vertices to which the external legs are attached; we then amputate both 
these legs and joint the vertices rj and p (the p vertex is always the first three- valent 
vertex in the rooted tree) by the propagator K(r),p). We took into account already 
the integration over £ thus obtaining an extra minus sign. We cannot perform this 
integration that easily in the second term in which the integration over £ is such that 
p > Cof. > V and the symbol J indicates that we must insert the integration 

/ d$ / dem', P )m,ri) 

Jp>C D ^>r) J oo a 

between the integrations over the variables p and rj. 

10 Conclusion 

We have defined quantum version of algebraic geometry notions, which allows us to 
solve the loop equations in the arbitrary /3-ensemble case. 
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The notion of branchpoints become "blurred", a branchpoint is no longer a point, 
but an asymptotic accumulation line along which we integrate instead of taking the 
residue at the branch point. 

Another surprising property pertains to the cohomology theory, which makes sense 
only if the cycle integral of any form depends only on the homology class of the cycle, 
i.e., we need all forms to have vanishing residues at the zeros Sj. This "no-monodromy" 
condition is automatically satisfied for our forms coming from the Schrodinger equation 
and it is equivalent to the set of Bethe ansatz equations satisfied by S{, similar to what 
takes place in the Gaudin model 

In contrast to paper I, there is no explicit dependence on the chosen sector; however, 
even the total number of ^.-cycles and rank of the period matrix may vary depend- 
ing on the choice of cuts in the complex plane. This might pertain to that we do 
not have actual finite-genus (classical) Riemann surface: analytical continuation may 
never result in sewing the corresponding solutions to the Schrodinger equation, and we 
therefore deal with different (finite-genus) sections of an ambient infinite-genus surface. 
So, indeed, the genus is no longer deterministic. 

Using the sectorwise approach, we can define the symplectic invariants; in Ap- 
pendix [D] below we present the first nontrivial calculation of this sort: the dependence 
of the leading term on the occupation numbers. 

In this paper, we restricted ourselves to the case of hyperelliptic curves, i.e. second 
order differential equations, that corresponds to 1-matrix model. The first straightfor- 
ward generalization pertains to including the logarithmic potentials into consideration, 
which would produce the Nekrasov functions nonperturbatively in the parameter 62/61. 
A more challenging problem is to generalize this approach to linear differential equa- 
tions of any order, which would correspond to a 2-matrix /3-ensemble model. In this 
case, we are also presumably able to define the notions of sheets, branchpoints, forms, 
and correlation functions. We expect also the preservation of the Bethe ansatz prop- 
erty ensuring a no-monodromy condition claiming that all cycle integrals depend only 
on the homology classes of cycles. The difference between the hyperelliptic case and 
the general case is comparable to the difference between the patterns of papers [9] and 
[6], i.e., the definition of the kernel K must be more involved and less explicit, but we 
postpone this discussion for further publications. 

It would be also interesting to see whether the quantities Th possess a symplectic 
invariance, or more precisely a "canonical invariance", i.e., whether they are invariant 
under any change (x, y) — > (x, y) such that [y, x\ = [y, x] = h. 
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A Appendix: Proof of theorem 15.2 

We now prove theorem 15. 2\ that all Wn s satisfy the loop equation, i.e., 



h i 



r=0 IC.J 

+ y dx ( W ( n h \x, J/{ Xj })5 a , aj - W ( n h) Cx v J/{ Xj }) \ 
j X] \ {x- •'', ) / 

is a polynomial in x of degree at most d — 2. 

From the definition, we have (with U from (15. 4p ) 

J) = ±£dzK@,z) (ui% 1 \z,z,J)+Y^B$,z)W^\z,^ . 

(A.2) 

Acting by H \ + d x j on K(x, z) gives l/{x— ^)+X)f=i hj( x )Cj(z), and the second 

part is obviously polynomial satisfying assertions of the theorem. Pulling the contour 
of integration w.r.t. z to infinity (with x originally outside the integration contour) 
and taking into account that the integral at infinity vanishes thanks to the asymptotic 
conditions, we find that only the residue at z = x and the residue at z = Xj in the 
second term in the brackets give nonzero contributions; the result of integration reads 

3 3 

so taking into account (I5.3p . we obtain the assertion of the theorem. □ 
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r(0) 

ClllUl)' Ul vv 

Theorem B.l The three-point function is symmetric 



B Appendix: The symmetricity of W% 



proof: 

Introducing Y := —2hip'/ip, is by definition 

1 



VK 



dx K(x°q, x)B(x\,x)B(x 2 , x) 



— (b dx Kn G\ G' 9 
4«7rjr CD 1 2 

-!- / dx AT ((^i + YK[ + Y'K x ){hK' 2 ' + + y'K 2 )) 

i'" .A-,., 

— / dx K ( ft 2 if"A~" + hY(K[K" + AT" AT') + hY'(K'(K 2 + K"K X 

+y 2 A^ 2 + yy'{K\K' 2 + a^ 2 ) + y' 2 ^^ ) 



where we have introduced a shorthand notation K p = K(x p , x), G p = G(x p , x), all the 
derivatives are w.r.t. x, and we omit indices indicating the sectors. 

The combinations K Q K x K 2 f(x), where f(x) is sector-independent (/ = 
1, U, U', . . .), vanish upon integration w.r.t. x because each of Ki is also sector- 
independent w.r.t. x. We can then use the Ricatti equation Y- 2 = 2KY- + AU to 
replace Y t 2 by 2KY( and YjY( by HY", which gives 

W^\xo,x x ,x 2 ) 

— / dx K Q (hY(K' 1 K" + K"K') + HY'(K'(K 2 + K"K X ) 
4«vr J Cd ' 

+2hY'K[K 2 + HY"{K X K' 2 + K[K 2 ) + Y' 2 K X K 2 



-!- / dx AT (fty(^i^)' + KY'(K X K 2 )" + KY"(K X K 2 )' + F' 2 ^^ ) 
= ^~<f dx Y' 2 K K X K 2 + H(Y"K (K X K 2 )' - (YK )'K[K' 2 - {Y 1 K Q )' {K X K 2 )') 
= -±-1 dx Y ,2 K K X K 2 - h((YK )'K[K 2 + Y' K' (K X K 2 )') 

4~ / rfx Y' 2 K K X K 2 - hYK' K[K' 2 - hY'(K K[K 2 + K' Q K X K' 2 + K' K[K 2 ) 



(B.2) 

This expression is readily symmetric in xo,a;i,X2 as claimed in theorem 15.31 □ 



49 



C Appendix: Proof of theorem 15.3 

Theorem Q Each wi 9) is a symmetric function of all its arguments. 
proof: 

The special case of wf was proved in appendix [B] above. The symmetricity of the 
two-point correlation function was proved in Theorem 14.91 

For technical reason, it is easier to proceed with the proof for nonconnected corre- 
lation functions. We introduce two types of them: 

• the correlation function 

k 

partitions i =l 

{/!,..., 4} Of I 

that comprises partitions of only stable (2hj + nj — 2 > 0) type with Ij ^ 0; 

• the correlation function 

k 

wf\i)= £ flW^Hlj) (C2) 

partitions 3=1 

{/!,...,/*} Of/ 

that admits also two-point correlation functions in the sums, 2hj+rij—2 > 0, 
with /,• ^ 0; 

The symmetricity of all W« (/) with s + 2h' < n + 2h obviously implies the 
symmetricity of W^ 1 (/) 

It is obvious from the definition that W^_ x (xq, x\, . . . , x^) is symmetric in 
Xi, x 2 , ■ ■ ■ , x n , and therefore we need to show that (for n > 1): 

5^(?o,*i, /) - wSi(x\,xl, J) = 0, (C.3) 

where J = \X2, ■ ■ ■ , x n |. 

The proof is by recursion on — \ = 2h — 2 + n. 

Assume that all and with 2h' + k-2<2h + n are symmetric. We have: 

(x° ,x\, J) 
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We first consider the product of functions KBK in the second term: recalling that 
B(x\,x) = d x (d x — 2 ^^ )K(x\,x) and integrating by parts, we obtain 



— ® dx K(x° , x)B(x\, x)K(x, y) 

711 Jc Vx >y 

— (h dxK' x (x ,x)K' x (x 1 1 ,x)K(x,y) 

2m Jc Vx >y 

- I dxK x (Z,x)K(Z\,x)2^K(x,y) 

™ JC^>y W 

dx K(x°q, x)K' x (x\, x)K' x (x, y) 



2ni 
1 



2iri 

4 



Cv x >y 



+_L I dxK(Z,x)K(x\,x)2^-K' x (x,y). (C.5) 
2™ Jc Vx>y if>(x) 

The first and the last terms in the right-hand side are already symmetric w.r.t. the 
replacement xq -h- x\ and we disregard them. Integrating by parts in the third term in 
the right-hand side we obtain one more symmetric term with K(x°q, x)K(x\, x)K xx (x, y) 
(which we can disregard as well) plus the term with K x (x°o,x)K(x\,x)K x (x,y). Com- 
bining the result with the second term, we obtain 

/ dxK' x (xlx)K(x\,x) ( dx + 2^)K(x,y) 
2™ Je Vx>y V W / 

= — <f dxK^ ,x)K{x\,x)l — + y^h fi {x)Cp(y)), (C.6) 
Zlxl Jc^yy \x-y p J 

where the integrand is sector-independent w.r.t. the variable x, so only the residue at 
x = y (with the minus sign) contributes in the second term in (1C.4|) . which therefore 
becomes 

dy2K' y (Z,y)K(x\,y)wi h - 1 1 \y,y,J). (C.7) 
For the first term in ( 1C.4I) . we use the induction assumption writing it in the form 



1 f , 1 



2m hv x >y 2 ™ JC V 



dx^—(k dy K(x° ,x)K(x\,y)x 
x (2B(x, y) 2 S h=1>n=2 + 4B(x, y)W^ l 1 \x, y, J)' + W^~ 2) (x, x, y, y, J)')(C.8) 



where the prime indicates that no propagators of B(x,y) type enter the expression, 
and no singularity occurs in the corresponding terms upon interchanging the order of 
contour integration w.r.t. x and y. The last term is again obviously symmetric w.r.t. 
the replacement x -H- x 1 . 

The skew-symmetric part in the middle term is one-half of the residue coming from 
the double-pole —l/(x — y) 2 in the expression for B(x,y) (it comes again with the 
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minus sign due to the choice of contour ordering), so we obtain 

■^—j dy2K' y (xl,y)K(x\,y)wj l b y 1 1 \y,y, J)' + symmetric term, (C.9) 



which exactly cancel the term in (jC.7p except the only case g — 1, n = 2 in which we 
use that B(x,y) = —l/(x — y) 2 + W 2 (x, y) as x — > y, so 

2B(x, yf = 2{x - y)- A - 4(x - y)~ 2 W { 2 0) {x, y) + regular, 

the most singular first term results in the integrand K"'K that is sector-independent 
and vanishes, whereas the second term produces 

— I dy 2K'Jx Q , y)K(x\, y)wf\y, y) + symmetric term, 
2 ™ Jcv y 



which kills the last remaining possible term in the expression (|C.7p . The theorem is 
proved. □ 

D Appendix: Calculating in the Gaussian 

case 

In this appendix, we calculate the singular part of the third derivative of J-q and 
integrate the answer to obtain the singular part of J-q itself. Although we explicitly 
calculate only the Gaussian model case, we propose the answer for the general model 
free energy singular part. 

In the Gaussian model case with the potential V(x) = x 2 , we have four sectors of 
solutions with the asymptotic directions ±oo, ±ioo. As the basic solutions we take 
ip+{x) and ip-(x) that decrease at the corresponding imaginary infinities -Moo and 
—ioo. The real axis then plays the role of the ^4-cycle whereas the imaginary axis is 
the £>-cycle. 

We are interested in evaluating the singular part of the third-order derivative ^^ a - 
Let us first find out about the origin of this singular behavior. Apparently, local 
singularities at finite t appear when the solutions ip + and ■?/>_ coincide, which happens 
when ip± := ip n = H n (ix)e x2 ^ 2h , where H n are the Hermite polynomials, and 

h 2 d 2 x ip n {x) = x 2 i) n {x) + (2n + l)hj} n {x). 

From Corollary 16. II we have 



1 



dtl (2nif J B J B J B 



dz 1 dz 2 dz 3 WP(zi, z 2 ,z 3 ), (D.l) 
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and using that (z\, %, %) — § Cv ^K(z\,^)B{^2,C)B(zt,^) and that no singular- 
ities appear when integrating over Z2 and z% we obtain that each integral gives, by 

a 

Theorem I4.8[ just v (£) with a = ± and 
with the normalization constant C$ such that 



u o(0 = ^o^TTV / r±{p)dp (D.2) 

±ioo 



[«b(0 - = 1- 



Note that even in the case where ip + = ip- the functions t>o(0 and t>o(0 differ because 
of different lower limits of integrations, their difference is just C ^|y ip 2 (p)dp, 
and the normalization constant Cq at ip + = ip- = ip n is 

• (D - 3) 

The remaining integral w.r.t. Z\ in (ID. II) develops singularity when ip± — )■ ^ n because 
the function if^,^) develops a logarithmic cut on the £>-cycle, and using the explicit 
form ( 14. f p for the i^-kernel (in this simplest case, K = K), we obtain 

w = Vk^L ^iwL^ (a (D - 4) 

where the contour goes between ±00 and ^finfty encircling the point p. The 
singularity occurs when p (and, correspondingly, z) tends to — 200 for f/> + and to -Hoo 
for ■?/>_; this singular part comes from the residue at £ = p, and we obtain that the 
expression in ( 1D.4I) is 

/"=R°° If/" 2 ± 

zl dz T~TTT~\ \ d P^l(p) v o(P) + regular 

and using the explicit expressions (ID. 21) and ( ID. 31) for t> , we obtain 
sing. 

V 0^0 / 



± 







_J ±ioo 



Cj, (D.5) 



where the singularity occurs at the upper integration limit for z and p when if)_ — > 
ip n , and the term in the square brackets is nonsingular in this limit, so we can replace it 
by its limiting value, which cancel exactly the corresponding term in the normalization 
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constant Co (see (ID.3j) ; the integral w.r.t. z and p can be separated, and we obtain the 
final expression 



-i 2 



dx 



(D.6) 



We now calculate t when ifj + ,ip_ — >■ ?/v. Choosing ip-(x) = ip+(x) f*^ ^r^y and 

taking into account that the number of poles of solutions outside the ^4-cycle is n, we 
have 

" +0 ° ^ - . /" + °° 1 



—/in + ft 
—/in + ft 



oo 
+oo 



-hn + ft 



1 



ioo ^(|) ^ Jo 



(D.7) 



The first integral in the denominator diverges as ip + — > ip n and denoting this integral 
as A, we have 



-hn + ft- 



oo ip'i_(z) 



L 







0(A 



-2\ 



•ioo V+(0 



Comparing this expression with ( 1D.6[) . we obtain 

1 



sing 



n e Z +i0 , 



(D.8) 



(D.9) 



ft(n + t /ft) 2: 

that is, this derivative has double poles with the coefficient 1/ft at all points t = —hn, 
n = 0, 1, . . .. A function that exhibits such a behavior is obviously T- function, so we 
have that, up to an entire function, 

. (&? \ 1 



sing. 



h 



[Logr]"(t /ft), 



(D.10) 



and, in turn, 



sing. F ~h 2 [j log r] (to/ft). (D.ll) 

Turning to the asymptotic behavior of f log T (sc) at large positive x we observe that 
the leading term is \x 2 log x, which is exactly what we might expect from the matrix- 
model-like arguments: we must be able to apply semiclassical approximation at large 
positive t /h, and in this regime we have the leading asymptotic behavior of Gaussian 
matrix model, i.e., sing.J-o — f^o^S^o) modulo polynomial terms (of order not higher 
than two). 

We may therefore put forward the following conjecture. 

Conjecture D.l The singular part of J-q for any potential Vd+\{x) has the form 
h 2 J2t=i |Lf l°gr](ei/^) where e, are the occupation numbers on the cycles A4. 
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